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Context. The most massive binary system Eta Carinae has been recently established as a gamma-ray source by the AGILE and Fermi- 
LAT detectors. The high energy spectrum of this gamma-ray source is very intriguing. It shows two clear components and a lack of 
any evidence of variability with the orbital period of the binary system. 

Aims. We consider different scenarios for the acceleration of particles (both electrons and hadrons) and the production of the high 
energy radiation in the model of stellar wind collisions within the binary system Eta Carinae with the aim to explain the gamma-ray 
observations and predict the behaviour of the source at very high gamma-ray energies. 

Methods. The gamma-ray spectra calculated in terms of the specific models are compared with the observations of Eta Carinae, and 
the neutrino spectra produced in hadronic models are confronted with the atmospheric neutrino background and the sensitivity of 1 
km^ neutrino telescope. 

Results. We show that spectral features can be explained in terms of the stellar wind collision model between the winds of the 
companion stars in which the acceleration of particles occurs on both sides of the double shock structure. The shocks from the Eta 
Carinae star and the companion star can accelerate particles to different energies depending on the different conditions determined by 
the parameters of the stars. The lack of strong GeV gamma-ray variability with the period of the binary system can be also understood 
in terms of such a model. 

Conclusions. We predict that the gamma-ray emission features at energies above ~100 GeV will show significant variability (or its 
lack) depending on the acceleration and interaction scenario of particles accelerated within the binary system. For the hadronic models 
we predict the expected range of neutrino fluxes from the binary system Eta Carinae. This can be tested through observations with the 
large-scale neutrino telescopes, which will support or disprove the specific hadronic models. 
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1. Introduction 

Eta Carinae is the highest mass binary system detected up to 
now. It contains a star with a mass estimated in the range 
80 - 12OM0 (Davidson & Humphreys 1997, Hillier et al. 2001) 
up to 2OOM0 (Kashi & Soker 2010). The basic parameters of 
this star remain uncertain. Its luminosity is equal to 4.5 x 1O''L0, 
and the surface temperature to 2 x 10"^ K. The radius estimates 
of Eta Carinae can differ by a factor of a few, between 40 - 180 
Rq (Hillier et al. 2001). The star creates a very strong wind with 
a terminal velocity of ~ 500 - 700 km s"'. The mass loss rate is 
estimated to be in the range from 2.5 x 10""^ Mq yr"' (Pittard 
& Corcoran 2002) up to 10"^ Mq yr"' (Hillier et al. 2001). 
Moreover, this mass loss rate seems to change rapidly (Martin 
et al. 2010, Corcoran et al. 2010). Eta Carinae is suspected to be 
a binary system with a period of ~ 2022.7 + 1.3 days (Daminelli 
et al. 2008), characterised by high eccentricity e ~ 0.9 (Nielsen 
et al. 2007) and the semi-major axis of a = 16.64 AU (Hillier et 
al. 2001). The parameters of the companion star are even more 
uncertain. It is a WR type or an O type supergiant with a lu- 
minosity 9 X 10^ L0, and a surface temperature ~ 4 x 10'* K 
(Verner, Bruhweiler & Gull 2005, Mehler et al. 2010). The ra- 
dius of the companion star estimated for the above parameters 
is ~ 1.4 X 10'^ cm. The mass loss rate of the companion star is 
estimated to be M, 



comp 



10 Mo/yr and the wind velocity to 
be ~ 3000 km s^' (Pittard & Corcoran (2002). The timing of the 
periastron passage is defined as To - 702452819.8 (Damineli et 



al. 2008). The binary system is immersed in the massive nebula 
(A/neb ~ 12 M0, Smith et al. 2003), which likely originated in 
the huge outburst observed in this system in the year 1843. The 
radius of the nebula is estimated to be 3.4 x 10'^ cm (Smith et 
al. 1998). 

The GeV y-ray emission has been recently reported from 
the direction of Eta Carinae by the AGILE telescope (Tavani et 
al. 2009). The y-ray flux was steady during the period of over 
one year corresponding to the pre-periastron passage except for 
a single flare lasting for about two days. The analysis of the 
whole Fermi-LAT data shows the y-ray source to be consistent 
with the AGILE discovery (Abdo et al. 2010, Walter, Farnier 
& Leyder 2010). These observations, covering also the period 
of the periastron passage, suggest a steady nature of the source. 
The y-ray spectrum extends up to ~100 GeV, showing two dis- 
tinct components. The first one is consistent with a power law 
spectrum (spectral index 1.6 ± 0.2) and a cut-off at 1.6 GeV. 
The second one is described by a power law with spectral in- 
dex ~ 1 .9 without evidence of a cut-off. The non-thermal nature 
of the source is also confirmed by the hard X-ray observations 
with the Beppo-Sax (Viotti et al. 2004), the INTEGRAL (Leyder 
et al. 2008) and the SUZAKU (Sekiguchi et al. 2009) sateflites. 
This hard X-ray emission creates an additional flat component 
above the extended soft thermal X-rays which dominate below 
1.5 keV. The X-ray emission from Eta Carinae in the 22-100 
keV energy range is very hard (differential spectral index equal 
to s 1 + 0.4, Leyder et al. 2008). More recent observations with 
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INTEGRAL show that the hard X-ray component is well de- 
scribed by a power law spectrum (with photon index 1.8) with- 
out any strong variability at the periastron passage (Leyder et 
al. 2010). 

In general, non-thermal high-energy radiation from massive 
binary systems is interpreted as radiation from particles accel- 
erated at the stellar wind shock (e.g. Eichler & Usov 1993, 
Benagha & Romero 2003, Bednarek 2005, Reimer et al. 2006, 
Pittard & Dougherty 2006 or Walter, Farnier & Leyder 2010). 
Both, leptons and hadrons can contribute to the y-ray emission 
from the binary systems in the considered scenarios. Particles 
might also be accelerated in the large-scale blast wave pro- 
duced during the Eta Carinae explosion in the year 1843 (e.g. 
Ohm, Hinton & Domainko 2010). It is argued that accelerated 
hadrons cannot find enough targets for efficient y-ray produc- 
tion. However, accelerated electrons are able to produce GeV 
y-rays in the inverse-Compton process by scattering infrared ra- 
diation from the nebula. 

In this paper we analyse different scenarios for the high- 
energy production in the general wind shock model with an 
application to the most massive compact binary system Eta 
Carinae. In contrast to previous works, we consider the accel- 
eration of particles in two shocks with different properties ap- 
pearing from the side of the Eta Carinae star and the companion 
star. 



2. General scenario 

We consider a simple scenario for the high-energy produc- 
tion within the binary system of two massive stars similar to 
what was already discussed in the papers mentioned above. 
The massive stars produce very strong winds, which collide 
within the binary system and create a double shock structure 
(see Fig. 1). The location of the shock within the binary sys- 
tem, Rsh, can be estimated for the known value of the coefficient 
rj = Mcompi^'comp/(A^^ECi'Ec), which describes the pressure ratio 
of the winds from both stars. For the Eta Carinae binary sys- 
tem the wind pressure from the side of the Eta Carinae dom- 
inates. The value of t] is estimated to be ~ 0.2 by Pittard & 
Corcoran (2002) for the following parameters of Eta Carinae, 
vec -500-700 km s^^ and Mec ~ 2.5 x lO""* Mq yr^^ and 
the companion star, Mcomp ~ 10^^ Mg/yr and the wind veloc- 
ity on Vcomp ~ 3000 km s"^ . This means that the shock is bound 
around the companion star but it is closer to the surface of the 
Eta Carinae star than to the surface of the companion star con- 
sidering their stellar radii. In the example modelling we present, 
we apply the above parameters for the stars and their stellar radii 
equal to /?EC = 1.2x10" cm and /?comp = 1.4x10^^ cm. The dis- 
tance of the shock from the companion star can be obtained from 
^comp ^ ^ ^ and from Eta Carinae ^D/(l+ yfq), 

where D is the separation of the stars. For the periastron passage, 
£) = (1 - e)fl * 1.7 AU. Then, * 6/?comp ~ 8.4 x 10'^ cm 

and R^^ « l-4/?Eta » 1.7 x lO'-' cm. For the apastron passage, 
the shock is located about an order of magnitude farther from 
both stars owing to the highly eccentric orbit. However, the frac- 
tions of their wind powers, which are intercepted by the shock 
structure, are comparable. The above mentioned estimates are 
obtained on the base of a simplified scenario. They may not cor- 
respond to the real situation as it occurrs in very close binary 
systems. At first, the stellar winds accelerate from the stellar 
surface. They are slower close to the stellar surfaces and reach 
the terminal speeds only at some distance from the stars. This 
may influence the localization of the shock within the binary 




Fig. 1. Schematic representation of the supermassive binary sys- 
tem Eta Carinae. Electrons (e) and hadrons (p) can be acceler- 
ated on both shocks to different maximum energies due to dif- 
ferent conditions at the shocks from the sides of the Eta Carinae 
star and the companion star. The details of the considered sce- 
nario are given in the main text of the paper. 



system. Secondly, the radii of the companion stars are not pre- 
ciously known. This means that in reality the distances of the 
acceleration region from the stellar surfaces may differ from the 
values estimated above. Thirdly, at the periastron passage the 
winds may not even balance. Owing to the instabilities in the 
wind, the shock surface can obtain a very complicated struc- 
ture or it may even collapse onto one of the stars (e.g. Parkin et 
al. 201 1). Therefore, the above estimated values should be taken 
with caution. They are applied as example values, but whether 
they correspond to the real situation in the Eta Carinae binary 
system is uncertain. 

It is believed that shocks produced in collisions of stellar 
winds are able to accelerate particles to relativistic energies. 
Note, however, that the conditions at the shocks from the sides 
of both stars can be quite different. The mechanical pressure of 
the winds has to balance at the shock, i.e., parts of the power pro- 
vided to the shock region from the sides of both stars are com- 
parable. However, other physical parameters such as the mag- 
netic field strength, the density of the wind or the wind veloc- 
ity, can differ significantly. Therefore, it is expected that spectra 
of leptons and hadrons accelerated at the shocks from the sides 
of both stars can have different properties because of the differ- 
ences in the efficiency of the acceleration process, magnetic field 
strengths and energy losses of particles that can determine their 
maximum energies. The density of stellar photons in both shocks 
is expected to be similar owing to relatively small thickness of 
the double shock structure. 

In principle both shocks can accelerate leptons and hadrons 
to different maximum energies. Particles accelerated in difi'erent 
shocks can produce radiation at difl'erent energy ranges and at 
diff'erent dominant radiation processes. We also expect that de- 
pending on the shock (from the Eta Carinae or the companion 
star), particles can already lose energy close to the shock (i.e. 
within the binary system) or escape to the surrounding large- 
scale nebula. Therefore, high-energy radiation with very diflFer- 
ent properties can be expected in this relatively simple scenario 
of the two stellar wind interactions. 
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3. Acceleration of particles 

In general, the acceleration rate of particles at the shock can be 
parametrised by 

/'acc = ^c£//?L«0.1^-5B GeVs-', (1) 

where f = 10"^f_5 is the acceleration parameter, E is the energy 
of particle (in GeV), /?l is the Larmor radius of particles in the 
magnetic field B (in Gauss), and c is the velocity of Ught. We 
can estimate the characteristic acceleration time scale of parti- 
cles from 

Tacc = £//'acc ~ 10£/(^-5B) S. (2) 

This acceleration process can be limited either by the advec- 
tion of particles along the shock surfaces with the stellar winds 
or by their energy losses through different radiation processes. 
The characteristic time scale for the advection process along the 
shock is, 

Tadv = 3/?sh/Vw ~ 3 X 10^/?13/V3 S, (3) 

where /?sh = 10'^/?i3 cm is the distance of the shock from the 
centre of the star and = lO'^vs cm s ' is the stellar wind 
velocity. The radiation energy losses of particles can determine 
their maximum energies. Below we estimate the energy loss time 
scales for electrons and hadrons and their maximum allowed en- 
ergies at the shocks from the sides of the Eta Carinae and the 
companion star. 



where Xo = 3.9 x 10 /pw cm is the electron radiation length in 
hydrogen, p„ = MI{AnR\v^) » 3 X 10'°M_4/?7|v;' cm"^ is the 
concentration of matter in the stellar wind, and fiie mass loss rate 
of the star is M = M© yr"'. 

The synchrotron and IC cooling time scales are clearly 
shorter than the advection time scale of electrons along the shock 
and their bremsstrahlung energy loss time scale for the applied 
parameters of the stars in the binary system. Therefore, the max- 
imum energies of accelerated electrons should be determined by 
balancing the energy gains from the acceleration process with 
the energy losses on the synchrotron and IC processes. The effi- 
ciency of these two processes depends on the energy density of 
the magnetic and radiation fields in the acceleration region. The 
synchrotron losses dominate over the IC losses in the T regime 
for the magnetic field £sh > m(TllRl)comp + iTt/R^hcV'^ G. 
For the part of the shock that is the closest to the stars at the pe- 
riastron passage, this limiting value is Z?sh > 160 G. Depending 
on the specific parameters of the companion stars, we expect a 
variety of radiation scenarios from the binary systems: 

1. the synchrotron energy losses always dominate over the IC 

process, 

2. the synchrotron energy losses dominate only at the high- 
est energies (whereas the IC scattering occurs in the Klein- 
Nishina regime), 

3. the IC energy losses dominate for the whole energy range of 
accelerated electrons. 



3.1. Electrons 

The synchrotron and ICS processes for the electrons are impor- 
tant. The cooling time scale of electrons on the synchrotron and 
IC (in the Thomson regime) processes can be estimated from, 
Tsyn/ic - E^lPsynjic, whcrc E^^ is the electron energy and Psyn/ic 
are the energy loss rates of electrons on both processes. These 
time scales are estimated as. 
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where ctt is the Thomson (T) cross section, Ub and t/rad is 
the energy density of the magnetic and radiation fields, me is 
the electron mass, /?sh is the distance of the shock from a spe- 
cific star expressed in units of the stellar radius of the specific 
star, and T = 10''r4K is the surface temperature of the specific 
star. The cooling time scale in the Klein-Nishina (KN) regime 
can be roughly estimated by introducing the energy of electrons 
into the above formula for the IC losses, corresponding to the 
transition between the Thomson and the Klein-Nishina regimes, 
eJ'™ = ml/i3kjiT) « 97 /T4 GeV. Then, the cooUng time scale 
of electrons in the KN regime in the radiation field of both stars 
is approximately given by 
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These radiation processes can occur in a different way on both 
sides of the double shock structure. 

By balancing the synchrotron energy losses with the acceler- 
ation energy gains, we obtain the limit of the maximum energies 
of accelerated electrons at the shock. 



s(6) 



The cooling time scale of electrons on the bremsstrahlung pro- 
cess is almost independent of the electron energy. It can be esti- 
mated from 



Tbr ~ Xo/C ^ 4.3 X 10*RI^V3/M.4 S, 



(7) 



l9Q(^-5/B)^'^ GeV, 



(8) 



On the other hand, by balancing the acceleration time scale with 
the IC energy loss time scale (in the T regime), we estimate the 
maximum energies of electrons as 



■ 4(^^5fish) 



1/2 



-1/2 



GeV. (9) 



To consider the high-energy processes in the Eta Carinae bi- 
nary system, we have to fix the basic parameter that describes the 
considered scenario, which is the strength of the surface mag- 
netic field of both stars. As an example, we use the value of 
2 X 10^ G for the surface magnetic field of the companion star 
and 200 G for the Eta Carinae star. Then, the magnetic field at 
the shock from the companion star is estimated to be B^J^™^ ~ 60 
G and at the shock from the Eta Carinae star as B^JJ"'' ~ 100 
G. We also fixed the values of the acceleration parameter ^ for 
both shocks by estimating them from ^ ~ (Vw/c)^. Because the 
velocity of the winds from both stars differs significantly, we 
obtain, ^ec ~ 5 x 10 and ^, 



comp 



10 . Because of large differ- 
ences between these acceleration coefficients and magnetic field 
strengths at the shocks, the maximum energies of the electrons 
accelerated on both shocks (from the sides of Eta Carinae and 
the companion star) differ significantly. 

The magnetic field strength in the winds of stars has a com- 
plicated dependence on the distance from the star. In a small 
region close to the stellar surface the magnetic field has a dipole 
structure, i.e., B(R) oc R-^ (USO & Mekose 1992). However, in 
most cases this region is very small and it can be neglected. At 
larger distances the magnetic field has a radial structure. Then, 
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its strength drops as B(R) oc R^^. This dependence dominates 
for distances characteristic for the periastron and apastron pas- 
sages of the stars in the Eta Carinae binary system. The above 
dependence of the magnetic field strength on the distance from 
the star has interesting consequences. The maximum energies 
of electrons, determined by the balance between their acceler- 
ation efficiency and the synchrotron energy losses, should in- 
crease proportionally with the distance of the shock from the 
stellar surface provided that the acceleration efficiency is inde- 
pendent ofthe magnetic field strength (see Eq. 8). Note however 
that at aertain distance from the stars the adiabatic losses can 
dominate over the hrotron energy losses. On the other hand, the 
mum energies of electrons, determined by the IC energy losses 
in the T regime, should stay independent ofthe distance from the 
star (see Eq. 9). We conclude that if the acceleration process is 
saturated by the synchrotron energy losses, the spectra of elec- 
trons should extend to higher energies for the parts of the shock 
farther from the star or for the apastron passage of stars when 
the whole shock structure is at a larger distance from the stars. 
However, the maximum energies should be independent of the 
shock location if the electron acceleration is saturated by the IC 
energy losses in the T regime. This conclusion is valid provided 
that the acceleration efficiency of electrons does not depend on 
the distance from the stars and the advection time scale of elec- 
trons is always longer than the energy loss time scales. 

3.2. Hadrons 

On the other hand, hadrons lose energy mainly on collisions with 
the matter of the stellar winds. The time scale for the energy 
losses on pion production in proton-proton collisions can be es- 
timated from 

Tpp = (o-ppA;cpwr^ ~6.3x 10Xi^'3/^-4 s, (10) 

where o-pp ^ 3 x 10"^*" cm^ is the cross section for p - p tt 
production, k = 0.5 is the in-elasticity coefficient in this colli- 
sion, and pw is the density of the stellar wind at the shock region 
defined above. This energy loss time scale for the shock from 
the side of the Eta Carinae star is estimated to be ~ 4 x 10"* s 
and from the side of the companion star as ~ 3.4 X 10^ s at the 
closest region of the shock to the stars and at their periastron 
passage. The time scales for energy losses by hadrons and their 
escape (advection) from the acceleration site do not depend on 
their energy. By comparing these time scales (Tpp = Tadv), we 
determine the condition under which hadrons can interact close 
to the shock within the binary system. The following condition 
has to be fulfilled 

M_4 > 0.2/?i3V^. (11) 

This condition is clearly fulfilled at the periastron passage of the 
binary stars in the region of the shock from the side of the Eta 
Carinae star but is not fulfilled at the shock from the side of the 
companion star. Therefore, we expect that hadrons accelerated 
at the shock from the side of the Eta Carinae star interact effi- 
ciently within the binary system but those accelerated from the 
side of the companion star escape from the binary system to the 
surrounding nebula and interact there with the matter of expand- 
ing winds and the matter expelled during the past outbursts in 
the binary system (i.e., in the years 1843 and 1890). 

The maximum energies of accelerated hadrons can be esti- 
mated by comparing the energy gains from the shock with en- 
ergy losses on hadronic interactions or the advection from the 



acceleration site. The comparison of Tacc with Tpp allows us to 
estimate the maximum energies of hadrons 

£p 6.3^_55/??3V3/M_4 TeV. (12) 

If the condition given by Eq. 11 is not fulfilled, then the maxi- 
mum energies of hadrons are determined by the escape along the 
shock. Then, the maximum energies are estimated as 

£p « 30^_5B/?i3/v3 TeV. (13) 

Hadrons are expected to be accelerated at the shock not only 
at the periastron passage but also at other phases of the binary 
system. The advection time of hadrons scales proportionally to 
the distance between the companion stars but the collision time 
of hadrons is only proportional to the square of the distance 
(owing to the dependence of density of the wind with the dis- 
tance from the star). We expect that hadrons are accelerated at 
the shock from the side of the Eta Carinae star to the maximum 
energies, which are independent on the phase of the binary sys- 
tem. However, their colUsion rate with the matter of the wind 
decreases at larger separation of the stars within the binary sys- 
tem. On the other hand, hadrons, accelerated at the shock from 
the side the the companion star and farther from the periastron 
passage (or at the parts of the shock laying father from the star), 
should be accelerated to lower maximum energies. This is a con- 
sequence of the mainly radial structure of the magnetic field 
around the stars within the binary system (B(r) oc R^^). As a 
result, hadrons injected into the Eta Carinae nebula will have a 
compUcated spectrum, which can only be limited by the extreme 
values. The more precise determination of this spectrum is diffi- 
cult due to the unknown dependence of the efficiency of hadron 
acceleration in different parts of the shock and at different phases 
of the binary system. Therefore, we calculate the expected radi- 
ation output below only from the hadrons at limiting distances 
from the star. 

4. Scenarios for gamma-ray production 

The compficated double shock structure of the colfision region 
of the stellar winds, the possibile acceleration of electrons and 
hadrons, and the possible domination of different types of radi- 
ation mechanisms, that are responsible for the saturation of the 
acceleration process of particles provide the reach range of sce- 
narios in which particles can produce high-energy radiation in 
the context of the massive binary systems. Here we consider the 
production of high-energy radiation in some possible scenarios 
in the context of massive binary system with an elongated orbit. 
They will be tested by the available high-energy observations 
of Eta Carinae. They also predict some specific features, which 
will allow us to distinguish between them with the future more 
systematic X-ray and GeV-TeV y-ray and neutrino observations. 

At first, we conclude that the hard X-ray emission detected 
by Beppo-SAX, INTEGRAL, and SUZAKU caimot be pro- 
duced by electrons in the synchrotron process. The maximum 
energies of synchrotron photons produced by electrons whose 
acceleration is saturated by the synchrotron energy losses can be 
estimated from 

e^me(£/ficr)re^ 1-3^-5 keV, (14) 

where ye is the Lorentz factor of electrons estimated from Eq. 8, 
and Bcr = 4.4 x 10^^ G is the critical magnetic field. Because 

the maximum energies of observed hard X-ray emission is at 
least s ~ 0.1 MeV, the acceleration efficiency has to fulfil the 
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condition > 76 (see Eq. HM . This condition cannot be ful- 
filled in the shocks either from the side from the companion star 
or from the side of the Eta Carinae star If the acceleration of 
primary electrons is saturated by the synchrotron energy losses, 
we should observe a strong synchrotron component in the Eta 
Carinae spectrum from the central source peaking at ~ 1 keV. 
This is in-consistent with the observations of Eta Carinae be- 
cause the detected soft X-ray emission has a thermal nature and 
is clearly extended. Therefore, we conclude that the hard X-ray 
emission cannot be produced by electrons in the synchrotron 
process, which are directly accelerated at the shocks. 

Below we consider the two most likely scenarios for the 
high-energy emission from the Eta Carinae binary system, i.e., 
the production of the whole y-ray spectrum only by electrons 
(our model A) and the production of the y-ray spectrum in a 
composite model, part of the spectrum by electrons and part by 
hadrons (model B). 

4.1. Model A: electrons 

The non-thermal multiwavelength spectrum measured from the 
Eta Carinae binary system has its maximum in the y-ray energy 
range. We expect that the energy losses of accelerated electrons 
in the synchrotron process are less than their energy losses in the 
IC process. As we noted above, both shocks present within the 
binary system are characterised by different values of the accel- 
eration efficiency (^comp ~ 10""^and^EC ~ 5 x 10"^) and different 
values of the magnetic field strength at the shock. We fixed the 
values for the surface magnetic field of the companion star as 
ficomp = 2 X lO'' G and of the Eta Carinae star as Z?ec = 200 G, 
and assumed a mainly radial structure of the magnetic field in 
the stellar winds within the binary system (i.e. B(r) oc r"^). With 
these assumptions, we estimated the magnetic field strength at 
the shock from the companion star to be B x 60 G, and from 
the Eta Carinae star as Z? « 100 G at the periastron passage. 
We compare the characteristic time scales important for accel- 
eration of electrons and their energy losses for the parameters 
mentioned above (see Fig. 2). The acceleration of those elec- 
trons at the parts of the shock that are the closest to the compan- 
ion stars at their periastron passage is considered. We find that 
electrons can be accelerated at the shock from the side of the 
companion star to the maximum energies E™" ^ 80 GeV. These 
energies are owing to the saturation by the synchrotron energy 
losses (Fig. 2). The saturation of electron acceleration by the 
synchrotron energy losses happens for energies close to the tran- 
sition between the dominance of the synchrotron and the IC en- 
ergy losses in the KN regime. Therefore, in spite of the saturation 
of electron acceleration by synchrotron process, these electrons 
lose energy mainly on the IC scattering in the KN regime (see 
Fig. 2). The acceleration process of electrons at the parts of the 
shock at larger distances from the stars should also be saturated 
by the synchrotron energy losses because their acceleration time 
scale depends on the distance from the stars in a similar way as 
their energy loss time in the KN regime (i.e. Tacc tic/kn R^, 
see Eqs. 2 and 6). However, the maximum energies of electrons 
accelerated at the parts of the shock at larger distances from the 
stars should increase proportionally to R (see Eq. 8). 

The maximum energies of electrons at the shock from the 
side of the Eta Carinae star are determined by the balance be- 
tween the acceleration process and the IC energy loss process 
in the T regime. For the parameters mentioned above, they are 
estimated to be /s"™ « 7 GeV (see Eq. 9). These maximum en- 
ergies should not change in the parts of the shock farther from 
the stars or at other phases because they are independent of the 




log(E / GeV) 

Fig. 2. Characteristic time scales for electrons in the accelera- 
tion process in the shock from the side of the companion star 
(dot-dashed line) and the side of the Eta Carinae star (solid), 
the synchrotron cooling time of electrons at the shock from the 
side of the companion star (dotted), the IC cooling time scale in 
the Thomson regime (dot-dot-dashed) and in the Klein-Nishina 
regime (dashed). The intersection of the acceleration time scale 
with the IC cooling time scale in the Thomson regime gives the 
maximum energies of electrons accelerated at the shock from 
the side of the Eta Carinae star. The intersection of the accelera- 
tion time scale with the synchrotron cooling time scale gives the 
maximum energies of electrons accelerated at the shock from the 
side of the companion star. The parameters defining the acceler- 
ation scenario are B™"'' = 60 G, B^^ = 100 G, ^°"p = lO"'*, 
and^EC 



distance from the stars (see Eq. 9). Electrons accelerated at the 
shock from the side of the Eta Carinae star lose energy mainly on 
the IC scattering of stellar radiation in the Thomson regime. The 
energy density of the radiation at the shock region from the Eta 
Carinae star and the companion star are comparable. Therefore, 
y-rays produced in the IC process should have a quite isotropic 
distribution. 

In conclusion, we expect the presence of two components in 
the electron spectrum, which are accelerated by different shocks 
created in the collision of the winds from the sides of both stars. 
These electrons can cool efficiently in the radiation field of the 
companion stars. These two components in the electron spec- 
trum are responsible in our Model A for the two-component y- 
ray spectra observed from the Eta Carinae binary system. 

4.2. Model B: electrons and hadrons 

In the second model we assume that the sub-GeV y-rays are 
produced by electrons as in the case of model A. However, the 
second higher energy component in the y-ray spectrum is ex- 
pected to be produced by hadrons (see also Leyder et al. 2010). 
Acceleration of hadrons at the shock from the side of the Eta 
Carinae star is saturated by their energy losses on pion pro- 
duction. Their maximum energies are estimated from Eq. (12) 
on /i™^" X 250 TeV. These hadrons cool close to the accelera- 
tion site within the binary system, producing high-energy y-rays 
and neutrinos. On the other hand, acceleration of hadrons at the 
shock from the side of the companion star is limited by their con- 
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vection along the shock. Then, the maximum energies of hadrons 
in the part of the shock that is the closest to the stars at the peri- 
astron passage, are estimated to be /i™" a: 5 PeV (see Eq. 13). 
These hadrons escape to the nebula surrounding the Eta Carinae 
binary system and eventually interact with the matter within the 
nebula, producing high-energy gamma-rays and neutrinos. In 
principle, hadrons accelerated in both shocks can contribute to 
the observed hard y-ray emission. However, their relative contri- 
bution is difficult to estimate because of it depends on the details 
of the acceleration process, which are mainly unknown. On the 
other hand, hadrons accelerated at the shock from the side of 
the companion star accumulate within the nebula. They interact 
with the matter of expending winds of the stars and the matter 
entrained within the nebula during the past outbursts. 

Note that in principle the efficiencies of the particle acceler- 
ation in both shocks can differ as can the efficiencies of the ac- 
celeration of electrons and hadrons. Therefore, a different com- 
bination of processes can produce a variety of radiation outputs 
from these binary systems. 

5. Absorption of gamma-rays 

The TeV y-rays produced within the compact binary systems 
should be absorbed in a complicated way in the soft radiation 
of luminous stars (e.g. Bednarek 1997). These absorption effects 
should be taken into account for the case of the y-ray production 
within ~10 stellar radii from both stars. Therefore, we calculated 
the optical depths for y-rays injected close to the shock region in 
the radiation field of both stars as a function of their propagation 
angle and for their different injection distances from the stars. 
We applied the parameters of both stars as mentioned in the in- 
troduction. The optical depths are shown in Fig. |3] The angles 
for the Eta Carinae star and the companion star are counted out- 
wards from the specific star Clearly even in the most favourable 
case defined by the angle a - 90°, the y-rays with energies 
above a few tens of Ge V should be efficiently absorbed (the opti- 
cal depth above unity). These absorption effects should be taken 
into account when comparing the y-ray spectra produced close 
to the periastron passage of the stars with the observations of 
the Eta Carinae binary system. However, closer to the apastron 
passage, most of the shock structure lies at distances larger than 
~ 10 stellar radii. At these distances the optical depths drop be- 
low unity (see Fig. 3 for the dependence of the optical depth on 
the distance). Therefore, it is expected that absorption of y-rays 
produced close to the shock region will result in a clear mod- 
ulation of the y-ray signal from the Eta Carinae binary system 
already at sub-TeV energies. 

6. Comparison with observations of Eta Carinae 

We discuss the radiation output expected in terms of the model A 
and B (defined above) for different locations of the acceleration 
site from the companion stars, i.e. for different parts of the shock 
and different phases of the binary system. The calculated hard X- 
ray to y-ray spectra are compared with the recent observations 
of the Eta Carinae binary system (described in the introduction). 

6.1. Model A 

In this model we assumed that the GeV y-ray emission is pro- 
duced by electrons accelerated at the shock from the side of the 
Eta Carinae star and the multi-GeV hard y-ray emission is pro- 
duced by electrons accelerated at the shock from the side of the 



companion star. As an example, we assume that electrons are 
accelerated with the power law spectrum defined by the spec- 
tral index ainj = 2 up to the maximum energies as estimated in 
Sect. 3.1. 

For electrons accelerated at the shock from the side of the Eta 
Carinae star, the main energy loss mechanism is the IC scattering 
in the T regime. This process is efficient enough already close to 
the acceleration site because the cooling time scale of electrons 
on the IC process is clearly shorter than their advection time 
scale (see Eqs. 3 and 5). As a result, the equilibrium spectrum of 
electrons reaches the spectral index ofeq = 3 and the y-ray spec- 
tra from these electrons should have a spectral index close to 2. 
However, at low electron energies, the cooling time scale of the 
IC process starts to become comparable to the cooling time scale 
of the bremsstrahlung process (see Eqs. 5 and 7). The transition 
between which of these two radiation processes dominates oc- 
curs for electrons with energies of a few tens of MeV. Therefore, 
the photon spectrum from the IC process should steepen at ener- 
gies somewhere in the soft y-ray to hard X-ray energy range. The 
precise break energy in the IC spectrum is difficult to estimate 
because of the uncertain density of the matter in the shock region 
(it varies by a factor of ~100, see Pittard 2009). For low-energy 
electrons, for which the bremsstrahlung energy losses dominate, 
the equilibrium spectrum should be close to the injection spec- 
trum a^^™ = ffinj. Therefore, the IC spectrum produced by these 
electrons should steepen to the spectral index ~1.5. We compare 
the hard X-ray to y-ray IC spectrum expected from the electrons 
accelerated at the shock from the side of the Eta Carinae star 
with the observations of this binary system in Fig. 4 (see solid 
curve in the left panel), y-rays from the bremsstrahlung process 
are produced with energies comparable to energies of their par- 
ent electrons. The power in this emission only contains a small 
part of the whole power of accelerated electrons. Therefore, its 
contribution to the spectrum observed from Eta Carinae is on a 
much lower level. 

On the other hand, electrons accelerated at the shock from 
the side of the companion star can reach energies of the order 
of ~ 100 GeV. The electron energies can be even higher, pro- 
vided that their acceleration occurs in the part of the shock at 
larger distances from the stars (Sect. 3.1). Electrons also lose 
energy in the IC scattering process, which in this case occurs 
in the T and in the KN regimes (the synchrotron energy losses 
are neglected, see Fig. 2 for details of the derivation of maxi- 
mum energies of accelerated electrons). We show that this IC 
emission can be responsible for the the hard multi-GeV y-ray 
spectrum observed from the direction of the Eta Carinae binary 
system. For the injection spectrum of electrons with the spec- 
tral index ffinj = 2, the IC y-ray spectrum is also produced with 
the spectral index ay - 2, both in the T and KN regimes (see 
e.g. Blumenthal & Gould 1970). However, as we showed above, 
y-rays produced at the part of the shock relatively close to the 
companion stars are efficiently absorbed in the stellar radiation 
field. We take these absorption effects into account by applying 
the optical depths for the y-rays in the radiation fields of the Eta 
Carinae and companion stars and using simple absorption law, 
oc exp[-r(£'y)]. Because the inclination angle of the Eta Carinae 
binary system is estimated to be ~ 40°, we show as an exam- 
ple the y-ray spectra with the absorption effects in the radiation 
of the Eta Carinae star for the angle a = 50° (see dashed and 
dot-dashed curves in left Fig. 4). The unabsorbed y-ray spec- 
trum (produced at distances >107?ec from the Eta Carinae star) 
are shown by the dotted curves. We conclude that the considered 
model predicts a clear modulation of the y-ray signal at ener- 
gies above ~100 GeV with the period of the binary system. The 



W. Bednarek and J. Pabich: High-energy radiation from the massive binary system Eta Carinae 



7 




log(E/MeV) log(E/MeV) log(E/MeV) 

Fig. 3. Optical depths for y-rays in the radiation field of the Eta Carinae star (left panel) and its companion star (middle panel) 
as a function of energy of injected y-ray photons. Specific curves show the results for the propagation angle of the y-ray photon 
(measured from the direction defined by the injection place and the center of specific star) equal to: a = 0° (dot-dot-dashed curve, 
direction outwards of the specific star), 30° (dot-dashed), 60° (dotted), 90° (dashed), 120° (soUd). y-rays are injected at the distance 
of 1 ARec from the surface of Eta Carinae and at the distance of 6/?comp from the surface of the companion star. The dependence of 
the optical depths on the distance from the Eta Carinae star is shown in the right panel for the injection angle of a = 60° and the 
distances IARec (dotted), 2Rec (solid), 3Rec (dashed), 5Rec (dot-dashed), and IORec (dot-dot-dashed), respectively. 



y-ray spectrum produced at the periastron passage by electrons 
accelerated at the part of the shock that are the closest to the Eta 
Carinae star should cutoff just below ~ 100 GeV. 

6.2. Model B 

In this model the GeV y-ray peak is also produced by electrons 
accelerated at the shock from the side of the Eta Carinae star (as 
in model A). But the multi-GeV hard y-ray emission is proposed 
to be produced by hadrons that are also accelerated at the shock 
from the side of the Eta Carinae star (see also Walter et al. 2010). 
At the periastron passage the hadrons accelerated at this shock 
can reach energies ~ 250 TeV for the magnetic field strength 
at the shock from the side of the Eta Carinae star (see Eq. 12). 
These maximum energies are independent of the distance from 
the stars in the case of the radial structure of the magnetic field. 
The number of colUsions of relativistic hadrons with the matter 
of the Eta Carinae wind, already close to the acceleration site 
within the binary system, can be estimated as, 

^ Tadv/Tpp ^ 4.8M_4/(V^/?13). (15) 

For the parts of the shock close to the Eta Carinae star at the pe- 
riastron passage, the collision rate is estimated to be ~15, and at 
the distance of 10/?ec it is still ~1.5. Therefore, we conclude that 
these hadrons are efficiently cooled already close to the shock 
within the binary system through the large range of the binary 
phases. As a consequence, y-ray spectra from the decay of the 
neutral pions produced by hadrons accelerated at the Eta Carinae 
shock should not depend strongly on the phase of the binary 
system. We calculate the spectra of y-rays and secondary lep- 
tons from the decay of the pions produced in hadronic collisions 
assuming the power law spectrum of hadrons with the spectral 
index ap = 2. In Fig. 4, y-ray spectra from decay of n" are com- 
pared with the y-ray observations of the Eta Carinae binary sys- 
tem. 

The y-ray spectrum from decay of pions clearly extends 
through the TeV energy range. However, as in the previous 
model, the TeV y-rays produced within ~10 stellar radii should 



be efficiently absorbed in the thermal radiation of the Eta Carinae 
star (Sect. 5). We take this absorption effect into account by 
showing in Fig. 4 the unabsorbed y-ray spectra (thin dashed 
curves) and absorbed spectra for the parameters of the binary 
system as for model A. The maximum energies of the acceler- 
ated hadrons should not depend on the distance from the stars for 
the radial structure of the magnetic ffeld, provided that the accel- 
eration coefficient does not depend on the distance (see Eq. 12). 
Therefore, at locations close to apastron, the y-ray emission ex- 
tending up to a few TeV is in principle possible. We conclude 
that the Cherenkov telescopes should detect a clear modulation 
of the y-ray signal at energies above ~100 GeV with the orbital 
period of the binary system because of the different absorption 
conditions within the binary system. 

The secondary leptons, from decay of charged pions pro- 
duced in hadronic collisions, efficiently lose energy in the syn- 
chrotron process, reaching the equilibrium spectrum with the 
spectral index or^^ = 3. They produce synchrotron radiation with 
the spectral index a^*^ - 2. The spectral index of the hard X- 
ray emission from the Eta Carinae binary system is not precisely 
known (Viotti et al. 2004, Leyder et al. 2008, Leyder et al. 2010, 
Sekiguchi et al. 2009). However, it is reported to be not far from 
the above value. The synchrotron emission from secondary lep- 
tons should extend through the hard X-ray energy range. The 
spectrum of secondary leptons, from decay of charged pions pro- 
duced in hadronic collisions, extends up to 

^max ^ E^^/iSfx) ~ 1 TeV, (16) 

where fi is the multiplicity of pion production by protons given 
in Orth & Buflington (1976). Electrons with these energies 
(^e = wJeTe) are able to produce synchrotron radiation at the 
shock from the side of the Eta Carinae star with energies up to 

£sy„ ~ me(5/5cr)r^ ~ 5 MeV, (17) 

i.e., consistent with the observations of the hard X-ray emission 
from the Eta Carinae binary system (see Fig. 4). The level of this 
emission should also be comparable to the level of the hard y- 
ray component from Eta Carinae, which is indeed the case (see 



8 



W. Bednarek and J. Pabich: High-energy radiation from the massive binary system Eta Carinae 




Fig. 4. X-ray and y-ray spectra expected for the model A (left) and model B (right) at the periastron passage of the stars in the 
binary system Eta Carinae compared with the high-energy observations of this system. The measurements in the hard X-rays by the 
INTEGRAL satellite (Leyder et al. 2010) are marked by squares and those in the y-ray energy range by the Fermi-LAT detector 
(Walter et al. 2010, Abdo et al. 2010) are marked by triangles. The results of the calculations are shown for the model A (left 
panel): y-rays from electrons accelerated at the shock from the side of the Eta Carinae star (solid) and from electrons accelerated 
at the shock from the side of the companion star: without any absorption (thin dashed), with absorption for the y-ray propagating 
at the angle equal to 50° (towards the star) at distance from the star D = 1.4/?ec (dashed), and at the same angle but for the parts 
of the shock at the distance from the star D - IQR^c (dot-dashed). For Model B (right panel): y-rays from electrons, accelerated 
at the shock from the side of the Eta Carinae star, scattering stellar radiation in the T regime (solid), y-rays from decay of neutral 
pions produced by hadrons accelerated at the shock from the side of the Eta Carinae star: unabsorbed spectrum (thin dashed curve 
extending to lower energies), spectra absorbed in the stellar radiation for distances from the star: D = 1.4/?ec (thick dashed) and 
D = 10/?EC (dot-dashed). Secondary leptons from decay of charged pions (produced in hadronic collisions) lose energy on the 
synchrotron radiation (dot-dot-dashed curve). The y-ray spectrum produced by hadrons, which are accelerated at the shock from the 
side of the companion star and produce pions in the interactions within the Eta Carinae nebula, are shown by the thin dashed curve 
extending to higher energies. 



Abdo et al. 2010 and Walter et. al. 2010). Therefore, secondary 
leptons can contribute to (or even be responsible for) the hard 
X-ray emission from the Eta Carinae binary system. The spec- 
tral index of the IC component in the hard X-ray energy range 
produced by primary electrons accelerated at the shock from the 
side of the Eta Carinae star is predicted to be clearly flatter than 
2 (see Model A). Therefore, precise measurements in this energy 
range should indicate the origin of this lower energy emission, 
i.e., either it is produced by primary electrons in the IC process 
or by secondary leptons in the synchrotron process. 

In principle, the hard y-ray component could be also pro- 
duced by hadrons accelerated at the shock from the side of 
the companion star. However, these hadrons escape to the Eta 
Carinae nebula from the binary system without significant en- 
ergy losses close to the acceleration site. The Larmor radii of 
hadrons escaping to the nebula are always much smaller than the 
characteristic distance scale of the system (which is of the order 
of the distance from the binary system) for the magnetic field 
around the massive stars with the toroidal structure at a distance 
above ~ 10.Rec (USO & Melrose 1992). Hadrons are frozen in 
the stellar winds moving outwards from the binary system. We 
can distinguish two regions within the nebula in which hadrons 
interact with the matter, i.e., the inner region where the density 
of matter is dominated by the direct stellar winds and the outer 
region where the density of the wind drops below the average 
density of the Eta Carinae nebula created during the outburst 
about 170 years ago. We assume that at larger distances from 
the binary system the winds from the Eta Carinae star and the 
companion star can mix through the rotation of the binary sys- 
tem. The characteristic collision rate of hadrons can be estimated 



from 

"^"^ PP ^r-HTnebPnebCTppC. (18) 

where i?bin is the radius of the binary system assumed to be equal 
to ~ 10'^ cm, /?neb = 3 X lO'^ cm, Pw('") is the density of the 
wind from the binary system (see above), Vw is the velocity of 
the stellar wind, pneb is the average density of the nebula ~ 10"* 
cm~^, and Tneb « 5 X 10^ s is the age of the Eta Carinae nebula. 
For the applied parameters, N'^°^ can reach the values of a few, 
i.e., hadrons can completely cool within the nebula. In principle, 
hadrons injected into the expanding wind from the star could suf- 
fer from adiabatic energy losses. To check whether this process 
can significantly change the energies of the injected hadrons we 
estimate the distance. Rout, from the binary system up to which 
hadrons interact at least once with the ni^tter of the wind. The 
following equation has to be checked: C °"' ^'^''^'^^^ dr = 1. For 
the example parameters applied in our calculations we estimate 
that hadrons should mainly interact relatively close to the binary 
system, where the density of the wind is still high. Therefore, 
adiabatic energy losses of hadrons through the expansion of 
the wind can be safely neglected. We conclude that hadrons in- 
jected into the nebula surrounding the Eta Carinae binary system 
should interact efliciently, producing y-rays and neutrinos. This 
conclusion is in contrast to the expectations of Ohm et al. (2010), 
who predicted inefficient interaction of hadrons in the large-scale 
nebula with relatively low average density. However, we con- 
sider the inner, denser nebula with the radius of a few times 
smaller than applied in their case. This inner nebula is domi- 
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nated by the matter from the dense wind of the Eta Carinae star 
produced at the present time. 

The maximum energies of hadrons escaping from the binary 
system can be as high as ~5 PeV (see Eq. 13) for the part of 
the shock that is closest to the companion star at the periastron 
passage. These maximum energies drop with the distance of the 
shock from the companion star if the radial dependence of its 
magnetic field reaches values of ~250 TeV for the distance of 
20/?EC- As an example, we calculate the y-ray spectrum expected 
from hadrons within the nebula with the power law spectrum, 
applying the spectral index equal to 2 and the maximum energies 
of hadrons equal to 5 PeV (Fig. 4). If hadrons are accelerated 
efficiently at the shock from the side of the companion star, a 
y-ray emission extending up to multi-TeV energies is expected. 
This emission component should be steady, i.e., independent of 
the phase of the binary system because the absorption efl'ects in 
the stellar radiation fields within the Eta Carinae nebula are not 
important. 

The total energy content in relativistic hadrons, which are 
responsible for the observed y-ray luminosity, can be estimated 
from, P^^^ » 3Ly''Tpp, where Ly'^ is the y-ray luminosity ob- 
served from the Eta Carinae binary system above ~10 GeV, Tpp is 
the cooling time scale for hadronic collisions given by Eq. (10), 
and the factor of 3 is there because only about one third of the 
produced pion energy goes to y-rays. For Ly'- » 2 x 10^"* erg s~^ 
derived for the distance to Eta Carinae equal to 2 kpc (Abdo et 
al. 2010) and the density of matter characteristic at the border 
of the binary system, Pw(^bin = 10'"* cm) ~ 10** cm"-', we esti- 
mate the total energy content in relativistic hadrons within the 
iimer nebula to be P^g^ ~10'^^ erg. Because the energy losses of 
hadrons are very efficient in this model, hadrons reach a steady- 
state distribution in the inner nebula. Then, the power radiated in 
y-rays is comparable to the power converted from the accelera- 
tion mechanism to relativistic hadrons, Lp ^ ll^^ . We conclude 
that to provide the hard y-ray luminosity observed from the Eta 
Carinae binary system, hadrons have to take only ~10~^ of the 
power of the wind from the Eta Carinae star, which is estimated 
to be Lw = Mv^ /2 * 2 x lO" erg s"'. 

Hadrons might also be accelerated at the large-scale blast 
wave produced during the past explosion of the Eta Carinae star 
(e.g. Ohm, Hinton & Domainko 2010). This injection process 
can provide an additional source of high-energy hadrons within 
the nebula which was not considered in the present paper. 

7. Neutrinos from Eta Carinae 

Hadrons accelerated at both shocks should also produce neutri- 
nos in collisions with dense matter of the Eta Carinae wind and 
the matter entrained within the Eta Carinae nebula during past 
outbursts. In model B, we propose that hadrons accelerated at 
the shock from the side of the Eta Carinae star can interact ef- 
ficiently with the matter of the stellar wind, producing y-rays 
via decay of pions from hadronic interactions through the large 
range of the binary phases. The maximum energies of acceler- 
ated hadrons are estimated to be ~ 250 TeV. They should not de- 
pend strongly on the distance of the shock from the Eta Carinae 
star. We calculate the muon neutrino spectra produced by these 
hadrons by applying the scaling break model, which is suitable 
for hadrons with the above estimated energies (Wdowczyk & 
Wolfendale 1987). It is assumed that hadrons have a power law 
spectrum with the spectral index equal to 2. Such a flat spectrum 
is motivated by the very flat hard GeV y-ray emission observed 
from the Eta Carinae. We compare the neutrino spectra with the 
atmospheric neutrino background and the sensitivity of 1 km^ 
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Fig. 5. Muon neutrino spectra produced by hadrons which are 
accelerated at the shocks from the side of the companion star up 
to energies of 5 PeV (solid curve) and the side of the Eta Carinae 
star up to energies of 250 TeV (dashed). Hadrons from the side 
of the Eta Carinae star interact with the matter of dense wind 
close to the shock region within the binary system. Hadrons from 
the side of the companion star escape from the shock region 
into the nebula and interact with the matter of expanding stel- 
lar wind and the matter within the nebula expelled in the past. 
The power law spectrum of hadrons is applied with the spec- 
tral index equal to 2. The neutrino spectra are normalized to the 
hard y-ray spectrum observed from the Eta Carinae binary sys- 
tem assuming its hadronic origin via pion production. The range 
of the atmospheric neutrino background (vertical and horizontal) 
within 1° degree of the source is marked by dot-dashed curves 
(Lipari 1993). The 1 yr sensitivity of the km^ neutrino detector 
(expected for the IceCube detector) is marked by the thin dotted 
line. 



neutrino detector of the IceCube type (see Fig. 5). Unfortunately, 
the level of this neutrino emission is comparable to the level of 
the atmospheric neutrino background. Therefore, these neutrinos 
will be difficult to detect even with a 1 km^ neutrino detector. 

Hadrons accelerated at the shock from the side of the com- 
panion star can also produce neutrinos in colUsions with the mat- 
ter of the Eta Carinae nebula. We show that this scenario can be 
tested by the future neutrino observations of the Eta Carinae bi- 
nary system. By assuming that y-rays, also produced by these 
hadrons, are responsible for the observed hard GeV y-ray emis- 
sion from the Eta Carinae binary system, we estimate the ex- 
pected muon neutrino spectrum from hadrons accelerated at the 
shock from the side of the companion star. The neutrino sig- 
nal should result as a contribution from hadrons accelerated at 
the shock during the whole range of phases of the binary sys- 
tem. Their maximum energies at these different locations vary 
between 5 PeV (periastron) to 250 TeV (at distance of the shock 
equal to D = 20/?ec)- Therefore, the total neutrino flux from 
hadrons, accelerated at the shock from the side of the compan- 
ion star and interacting with the matter of the Eta Carinae neb- 
ula, should be contained within the range defined by the dashed 
and solid curves in Fig. 5. There is a chance that these neutri- 
nos would be detected by the 1 km^ detector because the pre- 
dicted spectrum is above the atmospheric neutrino background 
and above the sensitivity of the 1 km^ neutrino detector at ener- 
gies between a few TeV and a few tens of TeV. Non-detection 
of the muon neutrino signal in this energy range (or detection 
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on lower level) will give us important information on the accel- 
eration efficiency of hadrons at the shock from the companion 
star 

8. Discussion and conclusion 

We investigated the consequences of the acceleration of particles 
(electrons, hadrons) at the shocks appearing within the binary 
system of two massive stars as a result of stellar wind collisions. 
In contrast to previous studies, we considered the acceleration of 
particles in two shock structures which appear on both sides of 
the contact discontinuity from the sides of the Eta Carinae star 
and the companion star. These two shocks have different proper- 
ties owing to differences in the surface magnetic field strengths 
(Bec - 200 G and Bcomp = 2 x 10^ G) and acceleration efficien- 
cies depending on the stellar wind velocities (^ec = 5 x 10"^ 
and ^comp - 10""^). The possibility of acceleration of particles in 
both shocks, either electrons or hadrons or both, can turn into 
a variety of radiation scenarios in which the high-energy y-rays 
can be produced. As an example, we considered in a more de- 
tail two general scenarios for some specific values of the surface 
magnetic fields of the stars and acceleration efficiencies at the 
shocks created from the sides of the Eta Carinae star and the 
companion star In the first one (model A), all high-energy emis- 
sion comes from electrons accelerated in these two shocks. In 
the second model (model B), soft y-ray emission comes from 
electrons and the higher energy y-ray component is produced by 
hadrons. Both models can successfully describe the y-ray spec- 
tral features of the Eta Carinae (e.g., two-component spectrum, 
lack of strong variability at the multi-GeV energies). However, 
the leptonic model predicts a clear variability of the y-ray emis- 
sion above ~100 GeV with the period of the binary system. This 
variability is caused by the strong absorption of the produced 
y-rays in the stellar radiation fields of both stars. On the other 
hand, TeV y-ray emission produced by hadrons in terms of the 
model B can either show variability in the case of their acceler- 
ation at the shock from the side of the Eta Carinae star (y-rays 
from within the binary system) or not show variability if they 
are accelerated at the shock from the side of the companion star 
(y-rays produced within the larger scale nebula surrounding the 
binary system). 

We have considered only a part of possible radiation scenar- 
ios that could be expected in this type of a wind collision model. 
For other values of the surface magnetic field of the companion 
stars, locations of the shocks, or acceleration efficiencies, the ac- 
celeration of electrons and hadrons can be saturated at other en- 
ergies. Therefore, a variety of spectral behaviours might be ex- 
pected from other massive binaries. The general spectral features 
of the Eta Carinae could also be investigated with other param- 
eters of the stars. Therefore, the set of parameters applied in our 
example modelling may not be unique. 

It is fairly difficult to conclude on the type of accelerated and 
radiating particles in the Eta Carinae binary system. In general, 
the leptonic model (model A) does not predict significant y-ray 
emission above a few TeV at any phase of the binary system. But 
the hadronic model (model B) predicts y-ray emission at ener- 
gies up to ~10-100 TeV. Moreover, y-ray spectra produced by 
hadrons injected into the Eta Carinae nebula should be indepen- 
dent of the phase of the binary system. The most evident confir- 
mation of the hadronic origin of y-ray emission could come from 
the detection of a neutrino signal from this source. Therefore, 
we calculated the neutrino spectra expected for the model B and 
compared them with the sensitivity of the 1 km^ telescope. We 
showed that the neutrino spectra produced by hadrons within the 



nebula can be above the atmospheric neutrino background and 
within the sensitivity of the 1 km^ telescope. However, the exact 
neutrino event rate in these telescopes is difficult to predict be- 
cause the energies of the injected hadrons depend on the phase 
of the binary system. The dependence of the injection rate of 
particles on the phase of the binary system cannot be reliably 
predicted at present. Only future precise measurements of the y- 
ray flux from the Eta Carinae as a function of the phase of this 
binary system may allow us to put limits on the injection rate 
of hadrons at different phases and thus allow us to estimate the 
expected neutrino event rates in large-scale neutrino telescopes 
more precisely. 

Acknowledgements. We would like to thank the anonymous referee for valu- 
able comments. This work is supported by the Polish MNiSzW grant N N203 
390834 and NCBiR grant ERA-NET-ASPERA/Ol/lO. Project co-funded by the 
European Union under the European Social Fund, POKL "Cztowiek najlepsza 
inwestycja". 



References 

Abdo, A.A. et al. 2010 ApJ, 723, 649 

Bednarek, W. 1997 A&A 322, 523 

Bednarek, W. 2005 MNRAS 363, 46 

Benaglia, R, Romero, G.E. 2003 A&A 399, 1121 

Blumenthal, G.R., Gould, R.J. 1970 Rev.Mod.Phys. 42, 237 

Damineh, A., Hillier D.J., Corcoran, M.F. et al. 2008 MNRAS 384, 1649 

Davidson, K., Humphreys, R.M. 1997 ARA&A 35, 1 

Eichler D., USO, V. 1993 ApJ 402, 271 

Hillier D.J., Davidson, K., Ishibashi, K., Gull, T. 2001 ApJ 553, 837 

Kashi, A., Soker, N. 2010 ApJ, 723, 602 

Leyder J.-C., Walter R., Rauw, G. 2008 A&A 477, L29 

Leyder J.-C., Walter R., Rauw, G. 2010 A&A 524, 59 

Lipari, R 1993 APh 1, 195 

Martin, J.C., Davidson, K., Humphreys, R.M., Mehner, A. 2010 AJ 139, 2056 
Mehner A., Davidson, K., Ferland, G.J., Humphreys, R.M. 2010 ApJ 710, 729 
Nielsen, K.E., Corcoran, M.F., Gull, T.R., Hillier, D. J., Hamaguchi, K., Ivarsson, 

S., Lindler D. J. 2007 ApJ 660, 669 
Ohin, S., Hinton, J.A., Doinainko, W. 2010 ApJ 718, 1610 
Orth, CD., Buffington, A. 1976 ApJ 206, 312 

Parkin, E. R., Pittard, J. M., Corcoran, M. F, Hamaguchi, K. 201 1 ApJ 726, 105 

Pittard, J.M. 2009 MNRAS 396, 1743 

Pittard, J.M., Corcoran, M.F 2002 A&A 383, 636 

Pittard, J.M., Dougherty, S.M. 2006 MNRAS 372, 801 

Reimer, A., Pohl, M., Reimer, O. 2006 ApJ 644, 1118 

Sekiguchi, A., Tsujimoto, M., Kitamoto, S., Ishida, M., Hamaguchi, K., Mori, 

H., Tsuboi, Y. 2009 PASJ 61, 629 
Sinith, N., Gehrz, R.D., Krautter J. 1998 AJ 116, 1332 

Smith, N., Gehrz, R.D., Hinz, PM., Hoffmann, W.F., Hora, J.L., Mamajek, E.E., 

Meyer, M.R. 2003 AJ 125, 1458 
Tavani, M., Sabatini, S., Pian, E. et al. 2009 ApJ 698, L142 
USO, V.V., Melrose, D.B. 1992 ApJ 395, 575 
Vemer E., Bruhweiler, F & Gull, T. 2005 ApJ 624, 973 
Viotti, R.F, Antonelh, L.A., Rossi, C, Rebecchi, S. 2004 A&A 420, 527 
Walter R., Famier, C, Leyder J.-C. 2010 A&A, submitted (arXiv:1008.2533"i 
Wdowczyk, J., Wolfendale, A.W. 1987 J.Phys. G 13, 411 



